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We constrain the possible time variation of the Higgs vacuum expectation 
value (v) by recent results on the primordial A He abundance (lp). For that, 
we use an analytic approach which enables us to take important issues into 
consideration, that have been ignored by previous works, like the w-dependence 
of the relevant cross sections of deuterium production and photodisintegration, 
including the full Klein- Nishina cross section. Furthermore, we take a non- 
equilibrium Ansatz for the freeze-out concentration of neutrons and protons and 
incorporate the latest results on the neutron decay. Finally, we approximate 
the key-parameters of the primordial 4 He production (the mean lifetime of the 
free neutron and the binding energy of the deuteron) by terms of ^ (where 
vq denotes the present theoretical estimate). Eventually, we derive the relation 

Y P ~ 0.2479 - 5.54 - 0.808 f 3 ^*) and the most stringent limit on a 

possible time variation of v is given by: —5.4 • 10~ 4 < < 4.4 • 10 . 



1. Introduction 



The standard model (lGrimthslll987l ) is a remarkably successful description of fundamen- 
tal particle interactions. The theory contains parameters - such as particle masses - whose 
origins are still unknown and which cannot be predicted, but whose values are constrained 
through their interactions with the so called Higgs field. The Higgs field is assumed to have 
a non-zero value in the ground state of the universe - called its vacuum expectation value v 
- and elementary particles that interact with the Higgs field obtain a mass proportional to 
this fundamental constant of nature. 

Although the questi on wh e ther t he fundamental constants are in fact cons tant, has a long 



history of study (see lUzan 



2003; Ichikawa & Kawasaki 



2002 



20031 ) for a review), comparatively less interest (TYbo fc Scherrer 



Kujat fe Scherrerll2000l ; IScherrer fc Spergellll993l ; iDixit &: Sher 
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19881 ) has been directed towards a possible variation of v. 

A macroscopic probe to determine the allowed variation rang e is given by th e netw ork of 
nuclear interactions during the Big-Bang-Nucleosynthesis (see lEidelman et al.l (12004 ) for a 
review of the Standard Big-Bang-Nucleosynthesis Model), with its final abundance of A He. 
The relevant key-parameters are the freeze-out concentration of neutrons and protons, the 
so called deuterium bottleneck (the effective start of the primordial nucleosynthesis) and the 
neutron decay. 

The major difference between our contribution and previous studies is, that we further im- 
prove this key-parameters by using an analytic Ansatz, exclusively. This analytic approach 
enables us to take important issues into consideration, that have been ignored by previous 
works, like the f-dependence of the relevant cross sections of deuterium production and 
photodisintegration, including the full Klein-Nishina cross section. Furthermore, we take a 
non-equilibrium Ansatz for the freeze-out concentration of neutrons and protons and incor- 
porate the latest results on the neutron decay. 

Finally we approximate the mean lifetime of the free neutron and the binding energy 



primordial 4 He abundance f 


Fukueita & 


Kawasaki 


2006; 


Olive & Skillman 


2004; 


Coc et al. 


2004; 


Izotov et al. 


1999: 


Izotov & Thuan 


2004; 


Luridiana et al. 


20031). 



We briefly note, that constraints on the spacial variation of v required a measurement 
of helium abundance anisotropy or inhomogeneity versus the position in the sky and an 
inhomogeneous theoretical BBN model. The homogeneous formalism used throughout the 
paper thus assumes a spacial invariance of the Higgs vacuum expectation value. 



2. Calculations 

All relevant processes of SBBN took place at a very early epoch, when the energy density 
was dominated by radiation, leading to a time-temperature relation for a flat universe: 



, 90fi 3 c 5 1 

where c is the velocity of light, k the Boltzmann constant, G denotes the gravitational con- 
stant and h is the Planck constant divided by 2n. g* counts the total number of effectively 
massless (mc 2 <C kT) degrees of freedom, given by g* = [g^ + !<?/), in which represents the 
bosonic and gj the fermionic contributions at the relevant temperature. The non-relativistic 



speci es are neglected, since their energy density is exponentially smaller (IKolb fc Turner 



19901 ). 



- 3- 



At very high temperatures (T ^> 10 10 K), the neutrons and protons are kept in thermal 
and chemical equilibrium by the weak interactions 

(1) n + e+ ^p + u e , 

(2) n + v e T± p + e~ and 

(3) n <=^p + e~ + z/ e , 

until the temperature drops to a certain level, at which the inverse reactions become ineffi- 
cient. This so called "freeze-out" -temperature TV and time tf de note the start of the effective 
neutron beta decay and detailed calculations (jMukhanovl 120041 ) derive 



Q 



Q 



+ 0.25 



7T- 



0.18\/-g*(T f ), 



(2) 



where Q is the energy-difference of the neutron and proton rest masses. At Tf the effec- 
tively massless species in the cosmic plasma are neutrinos, antineutrinos, electrons, positrons 
and photons. For the case of three neutrino families (jr& = 2 and gf = 10, which gives 
<7*(T» = 10.75. 



Following iMukhanovl (120041 ) we take a non-equilibrium freeze-out ratio of neutron number 
density (n n ) to baryon number density (un): 



n N 



exp 



-5.42(f^(T»)-i JZ(x + |) 2 (1 + e-~*)dx 



2y 2 (l + cosh(l/y)) 



dy, 



(3) 



where y = kTf/ Q. From now on the decay of free neutrons via n — > p + e + v e , with 



a mean lifetime (ISerebrov et al.l 120051 ) r n = 878.5 sec, can no longer be refreshed. Thus, 

t—t-f 

whereas the neutron density decreases as n n (t) = n n (tf) -e~^~ , the proton density increases 
as n p (t) = n p (tf) + (n n (tf) — n n (t)) and we obtain 



n 



n„ 



n N 



1 lj. 



(Tf)e 



- 1. 



(4) 



The next important step is the start of nucleosynthesis t N , usually referred to as the "deu- 
terium bottleneck" . The delay between tf and t^ is caused by the very low efficiency of direct 
production of light elements by successive collisions of several free protons and neutrons to 
one nucleus. In fact, nucleosynthesis proceeds through sequences of two-body reactions with 
the deuteron d as the intermediate product, via p + n — > d + 7. Accordingly, the small bind- 
ing energy of the deuteron ~ 2.225 MeV presents a severe problem for nucleosynthesis, 
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since energetic photons of the background radiation continuously disrupt the newly formed 
deuterons, until the temperature drops to a certain level T N , when the deuteron production 
gets the upper hand over photodisintegration. Because the decaying neutrons can no longer 
be refreshed by weak interactions after tf, the interval between tf and ijv plays an essential 
role for the outcome of the primordial helium production. 

Hence, we have to calculate the rates of deuteron production T( np ^ d ^, deuteron photo- 
disintegration Y^ dl ^ np ) and the expansion rate of the universe T exp , to determine t N respec- 
tively T N , when 



(np— >d"/) 
(dj^np) 



> 1 + 



r 



exp 



(d-y—tnp) 



(5) 



T exp for a radiation-dominated, flat universe is given by ^ with t from Eq ([T]). 

The rates for production and photo-disintegration of deuteron are given by the product of 

the relevant number density, velocity and cross section (cr): 



(np— >d7) 



n, 



&(np— >d-y) 



V 



(^(np—td'y) 



r 



(d-y^np) 



T^-y C &(dy— ►rap) 



(l + ^T)) £ c 



(6) 



where r\ ~ 6.14 • 10 -10 is the baryon to photon ratio based on WMAP (IHuey et al.l 120041 ) 
and n* denotes the number density of photons which effectively disintegrate the deuteron. 
These photons have to supply enough energy and must not loose this energy in much more 
likely Compton scattering with electrons instead of deuterons. 
The number density of photons at a certain temperature T is given by 



where ( is the Riemann zeta function. The number density of these photons, supplying a 
minimum energy E dis 3> kT is 



n (l>E dls ) 



8tt 



2 e , (kT x 3 

E 2 e~^dE = 8tt ( — 



E 



dis 



kT 



+ 1+1 



whereas only a fraction 



(d'y — 'rip) 



C r )' T (e7^e7) 



will successfully disintegrate a deuteron leading to 



n 



(■y>E dis )- 



@[drf— >rap) 



(9) 



We take ov 



( e7 ^ e7 ) ^ aKN(Ej) where cfkn denotes the Klein-Nishina cross section ( IRybicki h Lightman 
19791 ) for electron photon scattering and the mean incident photon energy E^ is given by 

oo 



En 



57T 



n 



{l>E dis ) 



(h c y 



E 3 e~^dE ~ E dis + kT. 



(10) 



E di . 
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The inter action cross section of p hoto-disintegration a^^np) can be derived by the 
calculations of iRustgi &: Pandeyl (119891 ) . We use a least mean square approximation within 
the incident photon energy range of 2.3 to 3.6 MeV and obtain 



57T 



^(d'y— >np) 



poo 

/ E 2 e~& [-2162.3 + 8. 1208 • 10 15 (£ + l.5kT)]dE (11) 

(12) 



n (i>E Ms ) (he) 3 

-2162.3 + 8.1208 • !0 15 (E dis + 2.5kT) {/ib}. 



The cross section for neutron capture <T( np ^ 7 ) is related to <7(d 7 -mp) by fjChen Savage 
19991 ): 



3^ 



^(d-y^np) 2m N c 2 (E 1 — B a 



(13) 



where we take E 7 from Eq ffTO]) . Collecting all terms, inserting into Eq (J5J) and taking into 



account, that E<n s = Bd — %kT we obtain an equation for T A 




(3/cTjv) 3 ff ( „)e^ 1 



m 3 c 6 



®{d~{— >np) 



*d 3 

Gg*(Tiv)fe 3 (T (e7 ^ e7 )A:T i ve fcT iV J 



90c 



( d'y — >■ np) 



10-3 4 (B d -ifeT w )2 



(14) 



where 



^(T^) = ^(T / )exp 



7?. 



90ft 3 c 5 



(15) 



and ^f-(Tf) is the reciprocal freeze-out ratio, given by Eq ([3]). The factor 10 34 on the right 
hand side of Eq (fill is due to the unit fib that we use for all cross sections a. 

The neutrinos have decoupled from equilibrium at about one MeV (above the rest 
mass energy of an electron) and thus before the annihilation of electron positron pairs. 
Therefore the entropy due to this annihilation is transferred exclusively to the photons, i.e. 
g*{T N ) ~ 3.36. 



Once sufficient deuteron has been produced, all other reactions 



\d + p 
ld+ld 
\He +1 He 



2-f/e and 
\He + 2p 
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proceed with significantly higher binding energies and the nucleosynthesis is no longer con- 
strained by photo-disintegration. It ends when the thermal energy is insufficient to permit 
the energetically favored synthesis reaction - the fusion of deuteron. 

With the assumption that at all available neutrons (as well as the same number of 
protons) have been synthesized to A He, which is not further transformed into heavier nuclei, 
because elements with nucleon mass number A=5 to A=8 are insufficiently stable to function 
successfully as intermediate products for nucleosynthesis at the available densities, we can 
calculate Yp, the final 4 He abundance by weight: 

2 _ 2gj-CZ» 2^(7)) 



Y P 



t N~ t 

e t " 



exp 



90fr 3 c 5 
327r 3 fc 4 G 



(16) 



where ^(7» ~ 
determines 



0.15709 is given by Eq © 
respectively. 



T f ~ f 0.64794 is given by Eq © and Eq <g 



rnum 
P 



For comparison with the most recent numerical result (ICuoco et al.l 120041 ) Yp 
0.2483, which assumes a mean neutron lifetime r n = 885.7 [s], we obtain (Eq [TB|) : Yp(r n = 
885.7) = 0.2483. 

For comparison with the observation-based result (jCoc et al.l 120041 ) Yp* s = 0.2479, we take 
JSerebrov et al.ll2005h r n = 878.5 [s] and obtain (EqUBj): Y P {r n = 878.5) = 0.2479. Further- 



more, this calculation shows, how sensitive Yp depends on the mean lifetime of neutrons. 

Taking into account our simple approach for the start of nucleosynthesis (EqEJ), where 
we neglected the fact, that the deuterons are not only destroyed by photo-disintegration but 
also consumed by the fusion of light elements, the concordance with the numerical as well 
as the observation-based result is very encouraging. 

This analytic expression for Yp therefore provides our basis for finding the dependence 
of Yp and the possible deviation of v from its present value vq, in order to finally constrain 
— by recent results on the primordial A He abundance. 



Crucial for the result of primordial nucleosynthesis is the moment or the correspond- 
ing temperature Tn, at which the production rate gets the upper hand. T/v depends on 
the binding ener gy of the deuteron B d . He nce we take the line ar fit of B d versus that 
has be en used by lYoo fc Scherrerl (120031 ) and iMiiller et al.l (120041 ) , based on iBeane fc Savage 
J2003h : 



B d (v) ~ B d (v ) (11 - 10 — 



m 



(17) 



where m n is the pion mass (the index again denotes the present value). As emphasized 
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by lYoo fc Scherrerl (120031 ). in our narrow range of interest, m\ oc v, leading to the final 
expression 

B d (v) ~ B d (v ) (11 - Wyfifa). (18) 

As Bd changes, Edi s , E 1 and the cross sections a^^np), ^(np^d-y) and cr( e7 ^ e7 ) change, 
accordingly. Concerning o"( e7 ^ e7 ) we furthermore have to consider, that the mass of the 
electron varies proportionally 

v 

m e (v) = m e (v ) — (19) 

which enters the Klein-Nishina cross section. 

Any variation of v, of course, changes the value of Q as well, but according to Eq (J2J), 
Tf is proportional to Q with the interesting consequence, that the freeze-out concentration 
of neutrons and protons does not change with a varying Q. 

By contrast, the Higgs vacuum ex pectation value defi nitely influences the mean lifetime of 



the free neutrons r n and following (jMuller et al.ll2004j ) we use the expression 



3.86 + 4 ^ + 1.52 me ~ me ° - 10.4^-7")-^™"°), (20) 

where a is the electromagnetic fine structure constant and and m u are the masses of the 
up- and down-quark (the index again denotes the present values). This approximation is 
the result of a linear analysis, based on the assumption, that only one single fundamental 
coupling changes with time while keeping the others fixed, respectively. Furthermore, Miiller 
et al. only consider standard model particles (with three neutrino families) to contribute to 
the energy density at BBN. 

Taking into account, that the elementary masses of the electrons and quarks linearly depend 
on v and disregarding the effect of a varying a (we take a as constant throughout this 
lettei0), we achieve 

r n (v) ~ r n (v ) (1 - 4.88 V -^-). (21) 

^0 

Finally, we derive a relation between Yp and v, to constrain the permitted variation of the 
Higgs vacuum expectation value by the primordial 4 He abundance: 

Y P ~ 0.2479 - 5.54 ' - 0.808 • (22) 



1 We also take H, k, c, G and rj as constant throughout the letter. 
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The predominant effect is the variation of the binding energy of the deuteron, followed by 
the mean neutron lifetime, whereas the changing mass of the electron with its consequences 
on the Klein-Nishina cross section (and the Compton scattering respectively) is almost neg- 
ligible. The relative weights can be quantified by the linearisation 

B d -B do \ fr n -T no \ „ ( m P 



Y P ~ 0.2479 + 0.105 ( — =° ) + 0.058 ( _ ) - 0.006 ( — ^ ) . (23) 

B do 



Results 



Using the observation based results of ICoc et al.l (120041 ) we derive 



-5.4 • 10" 



< 



v 



< 4.4 • 10" 



We avoid the term " observational results" because all cited publications more or less consist 
of interpretation of observational 4 if e-abundance plus theoretical input and constraints by 
the cosmic microwave background. Especially the different interpretation as a result of the 
deficiently understood systematics lead to incompatible data. Therefore, we separately state 
all publications and their constraints on v in the following table: 



Authors 


Y P 


Permitted variation — 

vn 


Fukueita & Kawasaki (2006) 
Olive & Skillman (2004) 
Coc et al. (2004) 
Izotov et al. (1999) 
Izotov & Thuan (2004) 
Luridiana et al. (2003) 


0.250 ±0.004 
0.2491 ±0.0091 
0.2479 ± 0.0004 
0.2443 ±0.0015 
0.2421 ±0.0021 
0.2391 ±0.0020 


(-8.0 ±I;g) • 10- 3 

(-1.5 tin) ■ io- 3 

(-0.49 ±ji) • IO" 4 
(4.3 t\'l) ■ IO" 3 
(6.8 + H) ■ 10- 3 

(io.i til) ■ io- 3 



Taking Vq—v as constant, the primordial A He abundance can be used for another inter- 
esting subject. Th e Higgs vacuum expectation value can be calculated theoretically by 



(IDixit fc Sherlll988l ) 



2-1/4^-1/2 



246.22 [GeV], 



(24) 



where Gf ~ 1.166371 [GeV~ 2 ] is the Fermi coupling constant (jEidelman et al.ll2005l ). based 
on measurements of the muon mass and lifetime. The uncertainty of determining t> , espe- 
cially the contributions of higher order terms, can now be constrained: 
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Authors 


I'll • S~1 ~\ T 

Compatible Vq m GeV 


Fukueita & Kawasaki (2006) 
Olive & Skillman (2004) 
Coc et al. (2004) 
Izotov et al. (1999) 
Izotov & Thuan (2004) 
Luridiana et al. (2003) 


245.56 t\H 
245.84 j£|jg 
246.21 ±g;il 
247.28 t° .tl 
247.90 t° H 
248.72 toil 



We thank K. Butler and A. Jessner for valuable comments on the manuscript. 
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